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ABSTRACT 

We apply the modified acceleration law obtained from Einstein gravity coupled to a 
massive skew symmetric field, F^\ , to the problem of explaining X-ray galaxy cluster 
masses without exotic dark matter. Utilizing X-ray observations to fit the gas mass 
profile and temperature profile of the hot intracluster medium (ICM) with King "/?- 
models" , we show that the dynamical masses of the galaxy clusters resulting from our 
modified acceleration law fit the cluster gas masses for our sample of 106 clusters with- 
out the need of introducing a non-baryonic dark matter component. We are further 
able to show for our sample of 106 clusters that the distribution of gas in the ICM 
as a function of radial distance is well fit by the dynamical mass distribution arising 
from our modified acceleration law without any additional dark matter component. 
In previous work, we applied this theory to galaxy rotation curves and demonstrated 
good fits to our sample of 101 LSB, HSB and dwarf galaxies including 58 galax- 
ies that were fit photometrically with the single parameter (M / L) sta rs- The results 
there were qualitatively similar to those obtained using Milgrom's phcnomcnological 
MOND model, although the determined galaxy masses were quantitatively different 
and MOND does not show a return to Keplerian behavior at extragalactic distances. 
The results here are compared to those obtained using Milgrom's phenomenological 
MOND model which does not fit the X-ray galaxy cluster masses unless an auxiliary 
dark matter component is included. 

Key words: dark matter — galaxies: clusters: general — galaxies: kinematics and 
dynamics — gravitation — X-rays: galaxies: clusters 



1 INTRODUCTION 

The question of galaxy rotation curves and cluster masses 
has been known to require some form of energy density that 
makes its presence felt only by its gravitational effects since 
Zwicky (1933) analyzed the velocity dispersion for the Coma 
cluster. By 1959 there were eight galaxy rotation curves 
available from radio observations demonstrating a "flat rota- 
tion velocity". In Brownstein & Moffat (2005) - our study 
of galaxy rotation curves without non-baryonic dark mat- 
ter - we demonstrated good fits to our sample of 101 LSB, 
HSB and dwarf galaxies including 58 galaxies that were fit 
photometrically with the single parameter (M/L) star s (29 
B— band and 29 if— band). As a follow-up, we apply here 
the same framework to the question of X-ray galaxy cluster 
masses (Moffat 2005a,b; Brownstein & Moffat 2005). 

In this framework, we apply a generalization of 
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Einstein's general relativity (GR) based on a pseudo- 
Riemannian metric tensor and a skew symmetric rank 
three tensor field F^x, called metric-skew-tensor- 
gravity (MSTG). A renormalization group (RG) frame- 
work (Reuter & Weyer 2004a,b) for MSTG was developed 
to describe the running of the effective gravitational 
coupling constant G, and the effective coupling constant 
7 C that measures the strength of the coupling of the 
F^vx field to matter. A momentum cutoff identification 
k — k(x) associates the RG scales to points in spacetime. 
For the static, spherically symmetric solution, the RG flow 
equations allow a running with momentum k and proper 
length £(r) — 1/k for the effective Newtonian coupling 
constant G — G(r), the coupling constant 7 C = 7 c (r), and 
the effective mass of the skew field /i = /i(r) where r denotes 
the radial coordinate. The form of G(r) as a function of 
r, obtained from the modified Newtonian acceleration 
law, leads to agreement with solar system observations, 
terrestrial gravitational experiments and the binary pulsar 
PSR 1913+16 observations, while strong renormalization 
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effects in the infrared regime at large distances lead to 
fits to galaxy rotation curves and X-ray cluster masses. 
In Moffat (2005b), a scalar-tensor- vector gravity theory 
is formulated, which yields the same weak field modified 
acceleration law as in MSTG Moffat (2005a), and equations 
that describe an effective running of G, 7 C and p with space 
and time. 

A fit to the bulk properties (the gas mass) of 106 X- 
ray galaxy clusters is achieved without exotic dark matter 
by a best fit result using a nonlinear least-squares fitting 
routine including estimated errors for the one free parame- 
ter in the RG flow equations, Mq, which fully describes the 
large distance renormalization of Newton's constant, Goo- 
The running of Newton's constant, G(r), is fully constrained 
by choosing the remaining free parameter in the RG flow 
equations, ro, to take on a value proportional to the scale 
of the galaxy cluster (for small clusters with galactic scale 
ICM gas masses) or to take on a particular constant value 
for the majority of the regular size galaxy clusters (where 
the ICM gas mass dominates over the masses of the individ- 
ual galaxies). Thus, the individual cluster mass profiles as 
functions of radial position throughout the range of X-ray 
observations are parameter free predictions; and are consis- 
tent with X-ray observed gas mass profiles without exotic 
dark matter. 

The galaxy cluster's bulk properties fit and the spe- 
cific fit to the individual clusters are compared to those 
obtained using Milgrom's phenomenological MOND model 
which does not fit the X-ray galaxy cluster masses unless 
an auxiliary dark matter component is included (Sanders 
2003). The & White (1988) were able to account for the 
MOND discrepancy between the X-ray observationally de- 
termined gas mass and the dynamical mass of the Coma 
cluster by increasing the MOND acceleration by a factor of 
four greater than that used to fit the galaxy rotation curves. 
Aguirre et al. (2001) present evidence from the central 200 
kpc of three clusters implying that MOND is inconsistent 
with the observed temperature gradient which inflates the 
discrepancy in the MOND acceleration to a factor of ~ 10. 
More recently, Pointecouteau & Silk (2005) use X-ray data 
from the XMM-Newton satellite for eight clusters of vary- 
ing temperature and masses to place constraints on the use 
of the MOND phenomenology. Without treating the MOND 
acceleration as a free parameter as opposed to a universal 
constant, MOND predicts dynamical masses greatly in ex- 
cess of the X-ray observations - necessitating the ad hoc 
addition of dark matter to explain away the missing mass. 
We are able to show that there is no missing mass when 
applying the MSTG acceleration law to galaxy clusters. 



2 ISOTROPIC ISOTHERMAL MODEL 

Recent observations from the XMM-Newton satellite sug- 
gest that the intracluster medium (ICM) is very nearly 
isothermal inside the region defined by the X-ray emission 
with temperatures ranging from w 1-15 keV (or 10 7 - 2 x 10 s 
K) for different clusters (Arnaud et al. 2001b). The combi- 
nation of the observed density profile, n e (r), and the tem- 
perature profile, T(r), obtained from X-ray observations of 
the galaxy cluster leads to a pressure profile, P{r), which 
directly leads to a mass profile, M(r), by assuming the gas 



is in nearly hydrostatic equilibrium with the gravitational 
potential of the galaxy cluster. Within a few core radii, 
the distribution of gas within a galaxy cluster may be fit 
by a King "/3-model". The observed surface brightness of 
the X-ray cluster can be fit to a radial distribution pro- 
file (Chandrasekhar 1960; King 1966): 
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resulting in best fit parameters, j3 and r c . A deprojection of 
the /3-model of equation (1) assuming a nearly isothermal 
gas sphere then results in a physical gas density distribu- 
tion (Cavaliere & Fusco-Femiano 1976): 

p(r) = p , (2) 

where p(r) is the ICM mass density profile. Provided the 
number density, n, traces the actual mass, we may assume 
that n(r) oc p(r), which according to Reiprich (2001) is ex- 
plicitly 
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and rewrite equation (2) 
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For a spherical system in hydrostatic equilibrium, the struc- 
ture equation can be derived from the collisionless Boltz- 
mann equation 
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where $(r) is the gravitational potential for a point source, 
a r and ae.4, are mass- weighted velocity dispersions in the 
radial (r) and tangential (6, <f)) directions, respectively. For 
an isotropic system, 

<J r = C0,0. (6) 

The pressure profile, P(r), can be related to these quantities 
by 



P(r) = <r r p(r). 



Combining equations (5), (6) and (7), the result for the 
isotropic sphere is 
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dr dr 
For a gas sphere with temperature profile, T(r), the velocity 
dispersion becomes 
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where k is Boltzmann's constant, p w 0.609 is the mean 
atomic weight and m p is the proton mass. We may now 
substitute equations (7) and (9) into equation (8) to obtain 
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Performing the differentiation on the left hand side of equa- 
tion (8), we may solve for the gravitational acceleration: 



a(r) = -- 



dr 
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pm p r 



d\n(p(r)) dln(T(r)) 
dln(r) 



dln(r 



(11) 



For the isothermal isotropic gas sphere, the temperature 
derivative on the right-hand side of equation (11) vanishes 
and the remaining derivative can be evaluated using the (3- 
model of equation (2): 



a(r) 



3f3kT 
pim v 



(r 2 +r 2 ) 



(12) 



3 MASS PROFILES 

Given an isotropic density distribution, p(r), the associated 
mass profile is 



M(r) = An p(r')r' 2 dr' , 



(13) 



where M(r) is the total mass contained within a sphere of ra- 
dius r. For the f3- model of equation (2), we may approximate 
the integral of equation (13) for r>r c and (3 < 1 (Reiprich 
2001): 
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This result clearly diverges in the limit as r — > oo; but galaxy 
clusters are observed to have finite spatial extent. This al- 
lows an approximate determination of the total mass of the 
galaxy cluster by first solving equation (2) for the position, 
r out , at which the density, p(r out), drops to w 10~ 28 g/cm 3 , 
or 250 times the mean cosmological density of baryons: 
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gas ~ 3(1 - /3) 

The dynamical mass in Newton's theory of gravitation 
can be obtained as a function of radial position by replacing 
the gravitational acceleration with Newton's Law: 

G M(r) 
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where Go is Newton's "bare" gravitational constant so that 
equation (11) can be rewritten as 
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and the isothermal (3- model result of equation (12) can be 
rewritten as 
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(19) 



Similarly, the dynamical mass in MSTG can be ob- 
tained as a function of radial position by substituting 
the MSTG gravitational acceleration law (Moffat 2005a, b; 
Brownstein & Moffat 2005): 
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so that our result for the isothermal /3-model becomes 
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We can express this result as a scaled version of equation 
(18) or the isothermal case of equation (19): 
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where we have made explicit the form of the running of G(r) 
as in Moffat (2005a,b); Brownstein & Moffat (2005): 
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In the limit of large r, we show in Moffat (2005a, b); 
Brownstein & Moffat (2005) that 
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and the total mass of the cluster in MSTG can be computed 
by taking equations (22) to the same limit: 

Mp \ 
Mmstg I 

It is a simple matter to solve equations (22) and (25) ex- 
plicitly for Mmstg(^) and Mmstg, respectively, by squaring 
both sides and subsequently applying the canonical solution 
to the quadratic equation. 

The derivation in Milgrom's phenomenological MOND 
model (Milgrom 1983; Sanders & McGaugh 2002) follows 
the same procedure, but utilizes the MOND gravitational 
acceleration law, described by 
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where n(x) is a function that interpolates between the New- 
tonian regime, p(x) — 1, when x 3> 1 and the MOND 
regime, p(x) = x, when i < 1. The function and criti- 
cal acceleration normally used for galaxy and cluster fitting 
are, respectively, 

p(x) = -j^=, (27) 

aoMiigrom = 1.2 x 10" 8 cm/sl (28) 

Applying equation (26) to either equation (11) or the 
isothermal case of equation (12) yields the MOND dynami- 
cal mass in terms of the Newtonian dynamical mass of equa- 
tion (18) or the isothermal case of equation (19): 

M N (r) 



MMOND(r) = 



1 + I lOMilgrom/ -p. 



(29) 



It is a simple matter to solve equation (29) explicitly for 
Mmond (r) by squaring both sides and subsequently apply- 
ing the canonical solution to the quadratic equation. 
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Unlike the Newtonian and MSTG isothermal spheres 
whose densities fall off as f /r 2 in the limit of large r , Milgrom 
(1984) showed that the MOND isothermal spheres have den- 
sities that fall of as r~ a where a m 4 at large radii. In 
the Newtonian and MSTG cases, the behavior is Keplerian 
at large radii and so the isothermal spheres have masses 
which diverge with r (Chandrasekhar 1960). Thus accord- 
ing to equations (17) and (20) the acceleration at large radii 
for Newtonian and MSTG isothermal spheres goes as 1/r. 
Conversely, the MOND isothermal spheres have convergent 
masses regardless of the spatial extent Sanders (2003): 
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The MOND acceleration law regularizes the integral of equa- 
tion (13) whereas a cutoff is necessary in the Newtonian 
and MSTG case, corresponding to the finite spatial extent 
of galaxy clusters as shown in equation (16). However, as 
clearly shown in Sanders (2003), the X-ray surface bright- 
ness distribution of these MOND isothermal spheres pro- 
vides a poor representation of the observed surface bright- 
ness distribution. This distinction will become apparent in 
the section ahead where we study the MSTG and MOND 
mass profiles of the Coma cluster. 



4 RUNNING OF NEWTON'S CONSTANT 

We have adopted the compilation of Reiprich (2001); 
Reiprich & Bohringer (2002) as our sample. The relevant 
cluster properties are listed in Table 1 arranged as follows: 
Column (1) lists the cluster name truncated to 8 charac- 
ters. Column (2) is the observed X-ray temperature. Col- 
umn (3) lists the ICM central mass density, po, of the ICM 
in units of 10 _25 g/cm 3 . Column (4) is the (5 parameter. 
Column (5) is the core radius parameter, r c , in units of kpc 



assuming Ho = 0.71_ 



(Eidelman et al. 2004). Column 



(6) is the position, r ou t at which the density, p(r ou t), drops to 
~ 10~ 28 g/cm 3 , or 250 times the mean cosmological density 
of baryons. Column (7) through (9) lists the ICM gas mass, 
Mgas, the Newtonian dynamic mass, Ms, and the MSTG dy- 
namic mass, Mmstg, respectively, each integrated to r out . 
Column (10) lists the "convergent" MOND dynamic mass. 

In order to calculate the MSTG dynamic mass we first 
need to phenomenologically determine the running of the 
parameters, Mo, and ro - this describes the running of New- 
ton's constant at the scale of clusters according to equation 
(23). However, unlike the case of the galaxy rotation curves 
where a satisfactory fit to LSB and HSB galaxy data is ob- 
tained with the parameters (Brownstein & Moffat 2005) 



M = 9.60 x 10 11 M , r = 13.92 kpc, 



(31) 



we found that a better fit was attained by dropping the 
simplifying assumption that Mo is constant across clusters. 
In fact, we were able to account for all dwarf galaxies smaller 
than 12 kpc in Brownstein & Moffat (2005) by allowing Mo 
to scale down to: 



M = 2.40 x lO n M0, r = 6.96kpc 



(32) 



(i) Mo oc Af" as , where M gas is the total ICM gas mass 
integrated to r out , 

(ii) Mo = constant, 

(iii) Goo = constant. 

It is clear that case (ii) is the limit of case (i) taking n — > 
and also that case (iii) is the limit of case (i) taking n — * 1. 
By plotting Mn against M gaa and then applying equation 
(25), we were able to constrain the parameter, Mo, using 
a nonlinear least-squares fitting routine including estimated 
errors . The fits for each of the three above schemes are 
shown in Figure 1, with the results as follows: 

/ \0.39±0.09 

(i) Mo = (58.8 ± 4.0) x 1O 14 M (to^%J 

(ii) M = (59.4 ±4.3) x 1O 14 M , 

(iii) Goo = 7.84 ± 0.27. 

The quality of each of these prescriptions for Mo is seen in 
their respective plots of Figure 1 where the MSTG mass is 
plotted against the ICM gas mass. Clearly, the plot corre- 
sponding to case (i) is the best scheme since the slope of 
unity best describes the bulk properties of the clusters. In- 
deed, the overall least sum of squares best fit is the scheme 
corresponding to case (i) with the other two acting as lim- 
iting cases. Meanwhile, both the MOND "convergent mass" 
and the Newtonian dynamic mass show a discrepancy with 
the ICM gas mass (Sanders 2003). 

In the galaxy rotation curves, we scaled ro down from 
equation (31) by a factor of 2 to account for the Dwarf galax- 
ies (Brownstein & Moffat 2005). For all but the smallest 
galaxy clusters where the mass of the ICM dominates over 
the mass of the individual galaxies, scaling ro up from equa- 
tion (31) by a factor of 10 leads to satisfactory fits, whereas 
for the smallest galaxy clusters where the mass of the ICM 
is of the order of galactic masses - with r ou t ^ 650 kpc - 
scaling ro to r ou t/10 leads to satisfactory fits. Thus our pre- 
scription for the running of Newton's constant (RG flow) for 
galaxy clusters is fully constrained: 

0.39 

M = 58.8 x IO^Mq ' 



(33) 



The three schemes we attempted to prescribe for the scale 
variation of Mo for clusters are as follows: 



ro = r out/10, r out < 650 kpc, 
ro = 139.2 kpc, r ou t > 650 kpc. 

5 INDIVIDUAL CLUSTERS 

We now turn to the individual mass profiles for each clus- 
ter. Since the MSTG prediction for the mass profile has no 
free parameters upon enforcing equations (33), we need not 
perform any nonlinear least-squares fitting for the individ- 
ual mass profiles as shown in Figure 2. The Coma cluster is 
shown first (with errorbands); and the remaining 105 clus- 
ters follow as thumbnails (without errorbands for clarity). 

The /3-model is not always a good fit to the data, how- 
ever, and so the mass profiles shown in Figure 2 have lim- 
ited reliability and the results of Table 1 show errors ranging 
from a few percent to in some cases 100%. Ota & Mitsuda 
(2002) show that some clusters observed with ROSAT and 
ASCA are expressed better by the double /3-model with a 
second core; while Pointecouteau et al. (2005) show that 
clusters observed with XMM-Newton are well described by 
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an NFW-type profile (Navarro et al. 1997) showing strong 
agreement with CDM simulations. While these more sophis- 
ticated mass profile models have their advantages, the sim- 
plest isothermal isotropic /3-model based upon hydrostatic 
equilibrium has the fewest parameters and we have used it 
to show that it is possible and meaningful to fit the X-ray 
galaxy cluster data without the need for exotic dark matter. 

There are also controversial measurements concern- 
ing the presence of temperature gradients in galaxy clus- 
ters (Reiprich & Bohringer 2002). XMM-Newton is ex- 
pected to probe the isothermal properties of the galaxy clus- 
ters. Early results such as the Arnaud et al. (2001a) study 
of the XMM-Newton observation of the Coma galaxy clus- 
ter and the Arnaud et al. (2001b) study of the XMM- 
Newton radial temperature profiles support the isothermal 
assumption deviating only near the cluster center where the 
observed X-ray temperatures decrease. 

In those cases where the uncertainty in the /3-model pa- 
rameters are sufficiently small (and therefore the /3-model 
provides a suitable fit), the MSTG mass profile traces the 
ICM gas mass profile down to sufficiently small radii and 
only deviates where the assumption that the ICM gas mass is 
isothermal is no longer valid. At sufficiently small radii, the 
observed X-ray temperatures decrease and thus the MSTG 
prediction based on an isothermal sphere is too large, con- 
firmed by Figure 2. 



6 CONCLUSIONS 

A gravity theory consisting of a metric-skew-tensor action 
(MSTG) that leads to the modified acceleration law (Moffat 
2005a, b; Brownstein & Moffat 2005) can be fitted to a large 
class of X-ray galaxy cluster mass profiles. The same accel- 
eration law can also be applied to fit a large class of galaxy 
rotation curves as in Brownstein & Moffat (2005) - which 
achieved one parameter fits for the case of the photometric 
observations with a single parameter (M/L) s tar S ; and also 
parametric fits which required no simplifying assumption 
on the constancy of (M/L) throughout the galaxy. Whereas 
MOND provides a good phenomenological fit to galaxy ro- 
tation curves, only MSTG fits both galaxy rotation curves 
and the X-ray galaxy cluster mass profiles without the need 
of introducing a non-baryonic dark matter component. 

In addition, whereas MOND necessitates the asymp- 
totic "flat rotation" velocity out to infinity, MSTG only gives 
the appearance of a "flat rotation" velocity within the con- 
fines of the galaxy and returns to the familiar Newton-Kepler 
form at large distances. Using the Sloan Digital Sky Survey 
(SDSS), Prada et al. (2003) have studied the velocities of 
satellites orbiting isolated galaxies. They detected approxi- 
mately 3000 satellites, and they found that the line-of-sight 
velocity dispersion of satellites declines with distance to the 
primary. The velocity was observed to decline to a distance 
of ~ 350 kpc for the available data. This result contradicts 
the constant velocity prediction equation of MOND, but is 
consistent with the MSTG prediction. 

The newest results from XMM-Newton reveal that the 
mass profiles of X-ray clusters show a steady rise out to the 
limit of X-ray observations (Pointecouteau et al. 2005) and 
do not show the convergent behavior of MOND isothermal 
spheres. However, the XMM-Newton results are consistent 



with the MSTG isothermal sphere whose densities fall off as 
1/r 2 without unseen non-baryonic dark matter. Further re- 
sults from XMM-Newton may reveal that the temperature 
profiles of X-ray clusters are not everywhere isothermal; and 
observed gradients may add a greater understanding of the 
physics of X-ray clusters of galaxies and the law of gravita- 
tion that best describes these largest of virialized objects. 
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Figure 1. A plot of the observed masses of the galaxy clusters in our sample listed in Table 1 versus the inferred Newtonian dynamic 
masses of the galaxy clusters and the best fit MSTG dynamic masses renormalized by three different schemes across all galaxy clusters: 
(1) Mq oc Mg as , where M gaa is the total mass of the cluster (black long dashed curve), (2) Mo = constant (blue dotted curve), and (3) 
Goo = constant (green short dashed curve). Red +-points are individual clusters from j3— model data (Reiprich 2001; Rciprich & Bohringcr 
2002). The MSTG dynamic mass in these three schemes is plotted in the remaining three graphs. Scheme (1) leads to the smallest sum 
of squares and is the overall best fit to the data with Mo ~ 58.8Af^ 9 [10 14 Mq]. Scheme (2) is an acceptable fit of the bulk properties 
with Mo ~ 59.4 X 10 14 Mq. Scheme (3) is a marginal fit of the bulk properties with Goo ~ 7.84Go- Both dynamic masses, Mmstg an d 
Mn are computed for isotropic, isothermal spheres cutoff at a spatial extent corresponding to a density of 10 — 28 g/cm 3 . Error bars are 
omitted for clarity — the computed errors are provided in Table 1 for completeness. 
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Figure 2. Galaxy Cluster Mass Profiles: Plot of the radial mass profile for the 106 clusters in the sample of Table 1. In all cases, 
the horizontal axis is the radius in kpc and the vertical axis is mass in units of 1O 14 M0. The red long dashed curve is the ICM gas 
mass inferred from X-ray observations according to the compilation of Reiprich (2001); Reiprich &: Bohringer (2002); the short dashed 
blue curve is the Newtonian dynamic mass; the dashcd-dotted cyan curve is the MOND dynamic mass; and the solid black curve is the 
MSTG dynamic mass. The Newtonian, MOND and MSTG dynamic masses are calculated within the context of the /3— model isothermal, 
isotropic sphere. The uncertainty in T, /3, r c and Ho is shown for the COMA cluster as a representative case of a cluster with a well fit 
/3-model; but error bands are omitted in the subsequent thumbnail plots for clarity. 
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Table 1. X-ray Cluster Properties of the Complete Sample: 
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Table 1 - continued X-ray Cluster Properties of the Complete Sample: 
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Note — We have adopted the compilation of Reiprich (2001); Reiprich & Bohringer (2002) as our sample. 
Column (1) Galaxy cluster name (truncated to 8 characters) Column (6) radius where Pgas ~ 10 — 28 g/cm 3 
Column (2) X-ray temperature Column (7) ICM gas mass integrated to r ou t 

Column (3) ICM central mass density Column (8) Newtonian dynamic mass integrated to ^ ut 

Column (4) /3-model f3 parameter Column (9) MSTG dynamic mass integrated to r ou t 

/3-model core radius parameter 



Column (5) 



Column (10) convergent MOND dynamic mass 



